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I. INTRODUCTION 

The question of how powerful quantum computers re- 
ally are remains to be answered. The difficulty of this 
question is not particular to quantum computational 
complexity. Classical question of whether there exists 
a polynomial algorithm for nondeterministic polynomial 
(NP) problems is one of the greatest problems in math- 
ematics. The prevailing opinion is that no such al- 
gorithm exists. Proving this seems to be exceedingly 
hard. One actually has no real idea of how to attack the 
problem. Due to an unintuitive character of quantum 
theory similar question for quantum algorithms seems 
only to be harder. Recently quantum adiabatic algo- 
rithm has been suggestcd for which the initial numeri- 
cal simulations showed polynomial scaling of the aver- 
age running time 0, Q for NP-complete problem. There 
are no known classical polynomial algorithms for NP- 
complete problems and some plausible arguments hint 
that it seems unlikely that a construction of quantum 
polynomial algorithm is possible 0. Nevertheless, even 
if the worst case complexity of quantum adiabatic algo- 
rithm is exponcntial, they might still provide a speed up 
for the average case performance. There were many sub- 
sequent numerical studies of the scaling of running time 
of adiabatic algorithm for different NP-complete prob- 
lems, some indicating exponential 0, ||, some polyno- 
mial dependence [||. While there exist analytic results 
for certain adiabatic algorithms (e.g. for Grover's search 
algorithm) 01111113, theoretical understanding of adi- 
abatic algorithms for NP-complete problems is still lack- 
ing. An exception is an analytical asymptotic expression 
for the energy gap which decreases exponentially for a 
particular choice of an initial Hamiltonian In view 
of the conflicting numerical results and in particular due 
to relatively small problem sizes amenable to numerical 
calculation theoretical understanding is greatly desired. 

Recently random matrix theo ry ( RMT) has been used 
to analyze adiabatic algorithm [12. IÏ3| , even though it 
has been noted 01 that it is not clear whether RMT 
applies to the low energy states. If RMT description 
would turn out to be applicable, we could use it to pre- 



dict the behavior of adiabatic algorithm for large prob- 
lem instances. In the present paper we are going to study 
statistical properties of eigenstates of adiabatic algorithm 
with the special emphasis on the question whether ran- 
dom matrix theory is applicable. While numerical re- 
sults for small problems in Ref. 0] supported the usage 
of RMT, we are going to show that for larger problems 
the behavior is quite different for problems having a non- 
degenerate ground state. 

In Ref. 12] the failure probability of quantum adia- 
batic algorithm is analyzed assuming random matrix the- 
ory spectrum, taking into account a cascade of Landau- 
Zener type transitions. From the result obtained, an ex- 
ponential scaling of the running time is suggested. They 
also numerically studied the distribution of ground state 
gaps for small 3-satisfiability (3-SAT) problems with 
ri = 8 variables. The distribution obtained showed a level 
repulsion which would support the usage of RMT also for 
the ground state. As we will show, the distribution ob- 
tained for the small n studied is not yet an asymptotic 
one and the behavior for larger n is very different. On 
the other hand, in Ref. |l3| the possibility of a polyno- 
mial running time is predicted, based on analysis of two 
RMT models, both giving essentially equivalent results. 
The main contribution to the failure probability in two 
models studied comes from the transitions to the bulk of 
the spectrum. 

In this paper we will study quantum adiabatic algo- 
rithm for 3-SAT problems having exactly one solution 
and show that while the bulk of spectrum is indeed well 
described by RMT the ground and the first excited state 
are far from RMT. In particular, it will be shown that the 
distribution of gaps does not show any level repulsion for 
sufficiently large n. RMT theory is therefore of limited 
use in describing the dynamics of Standard quantum adi- 
abatic algorithm for single-solution 3-SAT instances. As 
single-solution instances are thought to be the hardest, 
our results are important for the worst-case performance. 
What happens in the average case, when the number of 
solutions for certain vàlues of parameters may be large, 
remains to be explored. 
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II. QUANTUM ADIABATIC ALGORITHM 

We will study quantum adiabatic algorithm for 3-SAT 
with the Standard linear interpolation between the initial 
Hamiltonian H(0) and the final H(í), 

H(t) = (l-t)H(0) + tH(l). (1) 

The eigenstates of H(t) will be denoted by \ipi (í)) with in- 
teger index i denoting energy ordering, e.g., \ipo(t)) is the 
ground state. The energy and time t are dimensionless. 
The initial Hamiltonian H(0) is the sum of single-qubit 
Hamiltonians on each qubit, 

ff(0)=f>®l, A i=l(l 1 (2) 

i=l ^ ' 

while the final Hamiltonian H(í) is a sum of m three- 
qubit projectors, one for each clause, 

A m 

H(l) = -J2\C i )(C i \. (3) 

i=l 

A three-qubit projector given by Cj projects on the sub- 
space of states that violate i-th clause. The Hamiltonian 
H(í) therefore simply counts the number of clauses vio- 
lated by a given computational state. Somewhat uncon- 
ventional prefactor 4/a, with a — m/n, in H(l) is cho- 
sen in order to have time-independent trace, tr[if(í)] = 
Nn/2, with N — 2 n being the dimension of the Hilbcrt 
space. Two relevant parameters for 3-SAT are the num- 
ber of variables n and the ratio of number of clauses and 
variables, a = m/n. We used randomly generated 3-SAT 
instances having exactly one solution, i.e., the so-called 
single-solution random 3-SAT. We will predominantly fo- 
cus on instances with a = 3 as it has been recently nu- 
merically demonstrated Q that the gap for such prob- 
lems is much smaller than for those around the phase 
transition point for random 3-SAT |l4j . In addition, 
single-solution random 3-SAT instances with small a also 
seem to be hard for classical algorithms 0] . Because pre- 
vious studies focused on single-solution instances around 
the phase transition point we will for comparison occa- 
sionally also show the results for single-solution random 
3-SAT instances with a — 5, i.e., approximately at the 
location of the transition point for small n. 

To give an impression of how the eigenspectrum of 
H(t) OJ looks like, we show in Fig. ^ an example for 
n = 10. Two enlargements are shown, one for the central 
part of the spectrum and the other one for the lowest 
energies. Already at first sight one can see that there 
is a qualitative difference between both energy regions. 
While in the central part one has many avoided crossings, 
typical for a level flow described by RMT, in the lowest 
levels only one avoided crossing is prominent. In Ref. 5] 
it has been numerically shown that the failure probability 
is perfectly described by the Landau-Zener formula, tak- 
ing into account this single avoided crossing. One would 




FIG. 1: Spectrum for 3-SAT with n — 10 variables and a = 
3. The ground state gap occurring at tm'm = 0.562 has the 
value A = 0.111. The top figure shows enlargement of a tiny 
part of the spectrum in the so-called RMT core, while the 
bottom one shows the lower part of the spectrum. In the 
top figure, showing about 20 times smaller energy scales as 
the bottom one, there are many avoided crossings typical for 
RMT level flow, while in the bottom one only one avoided 
crossing between the two lowest levels is dominating. 



expect the eigenlevel statistics in the central part of the 
spectrum to be well described by RMT theory, i.e., we 
will have level repulsion due to many avoided crossings 
while on the other hand the level repulsion is expected 
to be very weak for the lowest levels (if present at all). 

RMT is successfully used to describe spectral statis- 
tic of complex systems, e.g., those with chaotic classical 
limit 0| . As the parameters of the system are changed, 
changing the dynamics from integrable to chaotic, the 
nearest neighbor spacing distribution changes from Pois- 
sonian to Wigner's surmise as predicted by RMT. Somc- 
times the perturbing parameter can be scaled out of 
the system and the transition from Poisson to RMT 
level spacing occurs as one goes from low to high en- 
ergies in the spectrum. Such is the case, for instance, 
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for the hydrogen àtom in a uniform magnètic field, see 
e.g., Ref. [TtJ. One of Wigner's main motivations to in- 
troduce RMT has been to describe resonances of neu- 
tron scattering on nuclei, i.e., the eigenspectra of nu- 
clei. Within Standard RMT the Hamiltonian has in- 
dependent matrix elements between all levels. To de- 
scribe excitations of nuclei it is actually more natural 
to use the so-called 2-body random matrix ensembles 
having only 2-body random interactions |l8j |. for a re- 
cent review, see Ref. 0. In particular, for 2-body 
random ensembles the distribution of spacings between 
the ground and the first excited state is more similar 
to semi-Poisson distribution p s p(s) — 4sexp(— 2s) than 
to Wigner's surmise for Gaussian orthogonal ensemble 
(GOE) pw(s) = fsexp(-s 2 7r/4) 20]. Note that for the 
semi-Poisson distribution there are more small spacings 
than for the GOE result. One can also argue that for a 
typical physical system the low energy spectrum will be 
dominated by some quasi-excitations, e.g., expanding the 
potential around the minimum one gets phonon-like ex- 
citations. Low energy spectral fluctuations as given, for 
instance, by the distribution of nearest neighbor spacings 
are therefore expected to be closer to those for integrable 
systems than to chaotic ones, i.e., more Poisson-like. Fur- 
thermore, it has been recently shown [2TL I22I ] that the 
degree of entanglement of the ground state for the quan- 
tum adiabatic algorithm for exact cover problem is much 
smaller than the maximal possible (as it would be, for 
example, for random vectors). 

All these results suggest that the low energy properties 
of quantum adiabatic algorithm could significantly devi- 
ate from those given by RMT. The aim of this paper is 
to show that this is indeed the case. But first, let us look 
at the high energy part of the spectrum where we expect 
RMT to hold. 



III. BULK PROPERTIES 




FIG. 2: (Color online) Violation of the PPT criteria for eigen- 
states, i = 0, . . . , 2™ — 1, at difïerent times í. Color (gray) en- 
codes the average minimal eigenvalue of p T/L , [Eq. J1J] . Inside 
the black circle the average value is larger than 0.15; we call 
this region the RMT core. All is for the same 3-SAT instance 
shown in Fig.Q 



We first test entanglement properties of eigenstates. 
For two qubits the positive partial transposition (PPT) 
is a necessary and sufficient condition for a state to be 
separable [23L |24| . For i-th eigenstate \ipi(t)) a t timc 
t the reduced density matrix pjk of the j-th and k~ 
th qubit is obtained by tracing over all other qubits, 
Pjk = ti£[\ipi(t))(ipi(t)\] , where £ is a set of all qubits 
apart from j-th and fc-th. We will denote the matrix ob- 
tained by partial transposition with respect to one qubit 
by pj£. If the smallest eigenvalue A m ; n (i,jfc,í) of pj£ 
is negative, the state pjk is entangled, otherwise, it is 
separable. We will use the smallest eigenvalue A m i n to 
measure 2-qubit entanglement of eigenstates. To obtain 
a quantity independent of two qubits j and k we will also 
average it over all pairs of qubits, 
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A mi n(ï,í) = — ^ -T V X min (i,jk,t). (4) 

nin — 1 ^— ' 

v ; iïk 

The dependence of A m i n (i,í) on the eigenvalue index i 
and time is shown in Fig. [21 If the eigenvectors are well 
described by RMT A m i n (i, t) should be large and positive. 
The reason is that for a random vector tracing over many 
qubits will result in a reduced density matrix that is very 
similar to the completely mixed one which is in turn sepa- 
rable. Therefore, whilc random vector on n qubits almost 
certainly represents an entangled state for any bipartite 
cut, say n/2 + n/2 qubits, it almost certainly does not 
give an entangled state when tracing over many qubits, 
specifically for a 2 x 2 degree of freedom reduced den- 
sity matrix. Numerical simulation for random states (ex- 
pansion coefficients are random Gaussian numbers) and 
n = 10 qubits gives A m i n (random) = 0.21, which is close 
to the largest vàlues attained for eigenvectors in Fig. [21 
The value A m i n (random) saturates for large n, while it 
weakly increases with n for smaller n. The asymptotic 
value is determined by the so-called induced measure on 
the space of density matrices |25f . Quantities other than 
Amin can also be considered and calculated for random 
2x2 degree of freedom matrices [2(| . In Fig. [21 one can 
observe that the largest vàlues of A m i n (i,i) are obtained 
in a certain "circle" of times t and eigenstates i. In Fig. [21 
we mark with a dashed curve the region where this value 
exceeds 0.15. We are going to call this central portion 
of the spectrum a RMT core, because eigenspectrum in 
this region can be well described by RMT. Note that in 
contrast to the RMT core there is a weak 2-qubit entan- 
glement present (negative A m ; n ) in the lower part of the 
spectrum (e.g., ground state). 

RMT behavior of states in the RMT core is confirmed 
also by studying bipartite entanglement. We divide n 
qubits into two halves and calculate the reduced density 
matrix p n /2 of the first n/2 qubits. The von Neumann 
cntropy of this reduced density matrix then characterizes 
bipartite pure state entanglement. The dependence of 
the entropy S(i, t) of the i-th eigenstate at time t, 

S(i,t) = -tr[/9 n / 2 hi2/0„/2] , (5) 
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FIG. 3: (Color online) The dependence of the entropy S(i, t) 
© of the reduced density matrix for the fïrst half qubits for 
eigenstates (index i) and time í. The entropy for maximally 
mixed state would be n/2 = 5. The black circle has the same 
location and the size as in Fig. [5] and the 3-S AT instance is 
also the same. 

is shown in Fig. [3J One can see the same structure as for 
the PPT criteria in Fig. [2 with the highest entropy eigen- 
states occurring at the same place as the highest vàlues 
of A min (the same black circle in the two figures) . For low 
energy eigenstates we again have strong deviations from 
RMT, i.e., small vàlues of the entropy. Remember that 
high vàlues of S(i,t) indicate strong bipartite entangle- 
ment while large vàlues of A m i n (z, t) indicate an absence 
of 2-qubit entanglement. 
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FIG. 4: Level spacing distribution for states from the RMT 
core. 9523 central eigenlevels of a single 3-SAT instance with 
n — 14 and a = 3 are used. Dashed curve is Wigner's sur- 
mise for GOE. In the inset a cumulative distribution is show, 
confirming the agreement with GOE also at small spacings. 

As a final test of RMT properties of states in the 
RMT core we studied the nearest-neighbor level spac- 
ing statistics, the paradigmàtic signature of RMT. We 
diagonalized the Hamiltonian at time t — 0.5 for onc 
3-SAT instance with n = 14 variables and a = 3, ob- 
taining all N = 16384 eigenvalues. For the central RMT 
core we choose to take 9523 eigenenergies in the range 
5.69 < Ei < 8.33 (the whole spectrum lies between 1.97 



and 13.00). Unfolding has been done by fitting a cúbic 
polynomial to the cumulative density in the used energy 
interval. The level density in this region is almost con- 
stant (» 2.5 • 10~ 4 in our case) and independent of a 
particular 3-SAT instance. The resulting level spacing 
distribution is shown in Fig.0|together with the Wigner's 
surmise pw( s ) = s7r/2 exp (— s 2 7r/4) for the GOE ensem- 
ble. One can see nice agreement with RMT also for small 
spacings seen in the cumulative distribution shown in the 
inset of Fig. 0] 

The situation is quite different in the low energy part 
of the spectrum. If using 1724 energies in the range 
2.5 < Ei < 5 (we start with the 8-th lowest energy), 
doing again cúbic unfolding, we get the level spacing dis- 
tribution shown in Fig.[S| While p(s) might seem to be in 
accordance with the GOE at first sight, the behavior of 
the cumulative distribution, shown in the inset, reveals 
that there are too many small spacings as compared to 
the GOE result. The level repulsion is therefore weaker 
in the low energy region. 
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FIG. 5: Level spacing distribution for low energy states (1724 
low energy states are used, starting with the 8th level) . In the 
inset, one can clearly see that there are more small spacings 
than for the GOE ensemble. 

The properties of eigenstates in the central part of 
the spectrum (RMT core) are therefore well described 
by RMT while there are deviations for low energy states. 
As the ground state and the first excited state are crucial 
for working of quantum adiabatic algorithm we will in the 
next section concentrate exclusively on the properties of 
the ground state. 



IV. GROUND STATE GAP 

A necessary condition for quantum algorithm to offcr 
exponential advantage over the classical is that the quan- 
tum states are sufficiently entanglcd, meaning that the 
entanglement, as quantified, for instance, by the maxi- 
mal Schmidt number, grows exponentially with size. If 
this is not the case, one could efficiently simulate quan- 
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tum evolution on a classical computer [27j. To describe 
the degree of entanglement we looked at the eigenvalues 
Xj, j = 0, . . . , 2™/ 2 — 1, of the reduced density matrix for 
the first n/2 qubits. Square roots of this eigenvalues are 
Schmidt coefficients for the n/2 + n/2 partition, 

2 „/2_ 1 ^ 

= E VA7h A )®K B >, (6) 

where and \x^) are the corresponding eigenvectors 
on the first and second n/2 qubits, respectively. How fast 
the eigenvalues Xj dccrcasc with j will teli us the degree 
of entanglement and if one can use efncient methods to 
simulate the evolution of such states (2^| . In Fig. H3 we 
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FIG. 6: Eigenvalues Xj, j = 0, . . . , 2 n/2 - 1 of the reduced 
density matrix for the first n/2 qubits ijüj. Pluses show data 
for a high lying eigenstate (i.e., 512í/i eigenstate) at í m i n = 
0.56 (the location of the minimal gap), crosses are for the 
ground state at í m i n = 0.56 while stars are for the ground 
state at t = 0.2. All is for the same 3-SAT instance from 

Fig. EJ 

show Xj for three different eigenstates: for the ground 
state and one high energy eigenstate from the RMT core 
at the location of the minimal gap and for the ground 
state at smaller time. One immediately notices the dif- 
ference between the ground state and the state from the 
RMT core. In the later the eigenvalues are much larger 
and decrease with j very slowly. On the other hand, for 
the ground state the eigenvalues Xj decrease much faster. 
This fast decrease of Xj for the ground state makes it 
possible to simulate ground state dynamics (e.g., quan- 
tum adiabatic algorithm) on a much smaller space than 
the full N dimensional Hilbert space. This has been ex- 
ploited to perform numerical simulation of the quantum 
adiabatic algorithm for much larger n than possible with 
the conventional methods |2l|. From Fig. [ü]we can, for 
instance, see that if we are content with the precision of 
say 10~ 3 one has to take ~ 25 eigenvectors in the case 
of the "random" state from the RMT core, while on the 
other hand we would need only ~ 5 eigenvectors for the 
ground state. Note that Xj are simply connected with the 



entropy S [Eq. (0] shown in Fig. For the three cases 
shown in Fig. [ü]we get entropies S(i = 0,0.56) = 0.90, 
S(i = 0,0.20) = 0.15 and S(i = 512,0.56) = 4.08. The 
entropy of the ground state 5(0, í) attains its maximal 
value at the position of the minimal gap j^. It is sig- 
nificantly smaller than for the excited states from the 
RMT core but still grows linearly with n, preventing the 
efficient classical simulation of the quantum adiabatic al- 
gorithm for NP-complete problems [22^ . 

While the entanglement of the ground state determines 
how efficiently we can classically simulate such an algo- 
rithm, the running time of the quantum adiabatic algo- 
rithm is predominantly detcrmined by the minimal gap 
A between the ground and the first excited state. The 
necessary running time for the wanted precision at the 
end can be simply determined trough the Landau-Zener 
formula [ 3. 

In Ref. [ïí| the authors found a GOE-like distribution 
for the minimal gaps A of small 3-SAT instances with n = 
8 variables. If such behavior would persist for larger n 
this would be advantageous because due to level repulsion 
we would have fewer small spacings, i.e., the running 
time could be smaller. One should note that here we are 
talking about the distribution of ground state gaps A for 
different 3-SAT instances, i.e., the distribution is mcant 
over the single-solution random 3-SAT ensemble (many 
spectra), whereas in the previous section we looked at the 
distribution of spacings within a single spectrum. Here 
we also do not do any unfolding as there is no obvious 
unfolding procedure for the lowest state. 

n=8, a=3 n=8, a=5 
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FIG. 7: Histograms of the distribution of minimal gaps A for 
1000 single-solution 3-SAT instances for different n and a. 
Spacings s are expressed in terms of the average spacing A 
that can be read from Fig. |Ü| For small n there are few small 
spacings, whereas for larger n their number increases. This 
transition happens sooner for smaller a. 

In Fig. [7| we show the distribution of gaps for 1000 
random 3-SAT instances at two different n and a. The 
gaps are expressed in terms of the average gap, s = A/A. 
We can see that for small n we indeed have a sort of 
level repulsion. For larger n = 16 the character is quite 
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different though. For a — 3 the distribution is more 
Poisson-like with an abundancc of small spacings. For 
a = 5a similar behavior can be observed, but it seems 
that for larger a the change from GOE-like to Poisson- 
like distribution takes place at larger n. 

The histogram for the largest case of n — 18, a = 3, 
we generated is shown in Fig. We can see that the dis- 
tribution is close to Poissonian with two important dif- 
ferences. There are more small spacings and more large 
spacings than one would expect for an exponential distri- 
bution. Of course, the probability to have zero spacing 
is zero, so for very small spacings p(s) goes towards zero, 
e.g., for the case in Fig. [5] cumulative distribution (not 
shown) grows as ~ s 1 ' 5 for very small spacings of order 

s » ícr 3 - ícr 2 . 
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FIG. 8: Level spacing distribution for 1000 single-solution 
random 3-SAT instances with n = 18 and a = 3. The full 
line is the exponential curve (Poisson spacing distribution), 
the arrow shows the position of the median spacing, s mct j = 
0.56, and the vertical line next to the median is the value 
of the median for an exponential distribution. All spacings 
are expressed in terms of the average spacing, s — A/A. No 
RMT-like level repulsion is present. 

In Fig.|H|we show in the main plot the median spacing 
for a = 5 and a = 3. Instead of showing A directly we 
divided it with n because the trace of H(t) oc n and so we 
expect that the eigenvalues themsclves will grow propor- 
tionally to n. We can see that the dependence of A/n on 
n is exponential for a = 3, whilc the asymptotic behavior 
for a — 5 has possibly not yet been reached. Such expo- 
nential decrease of the gap, suggesting exponential run- 
ning time, has already been numerically found in Ref. 0- 
Interestingly, the asymptotic decay rate for a = 3 agrees 
with 2 A/n x í/^/N (shown with dashed line), which is 
the same as the analytical asymptotic result obtained in 
Ref. [Ül for the initial Hamiltonian being a projector to 
the ground state. It might be that the worst-case per- 
formance (i.e., for small oi) is A ~ 1/yN regardless of 
the choice of H . Still, this issue needs to be explorcd 
in more detail. In the inset to Fig.[ü]we show for a = 3 
also the minimal, maximal, and the average A/n. While 
all seem to have exponential dependence on n, their de- 



cay rate is different. This could hint that by increasing 
n one gets increasingly more small and more large spac- 
ings. The same conclusion has been reached from the 
distribution of gaps in Fig. [8] where there is also a dif- 
ference between the median and the average A/n. The 
distribution of gaps therefore seems to change character 
with n and one cannot claim that the distribution found 
for n — 18 (Fig. [SJl is already the final asymptotic distri- 
bution. In fact, there even might not exist any stationary 
asymptotic distribution. Whcther these múltiple scales 
of the gap distribution can be remedied by some unfold- 
ing procedure is not clear. In any case the distribution 
of gaps is far from RMT prediction and, unfavorably, we 
have many small spacings that will necessitate large run- 
ning times of quantum adiabatic algorithm. 
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FIG. 9: The dependence of the median gap A/n on n for 
a = 3 and a = 5. Exponential dependence can be seen for 
a = 3 while for a = 5 the large n dependence is hard to 
infer (fitting a power law gives A ~ n~ ' 6 ). Dashed line is 
A/n = 2 ^jy ■ In the inset we show for a — 3 in addition to 
the median also the average, minimal, and maximal gap, all 
out of 1000 instances. 



V. PROBABILITY FLOW 

So far we have identihed two regions in the spectrum of 
quantum adiabatic algorithm for 3-SAT. The bulk prop- 
erties in the RMT core are described by RMT while the 
ground state properties and the ground state gap A are 
far from RMT predictions. In the present section we are 
going to explore how important the RMT core states are 
for the success of adiabatic algorithm. 

Let us denote the solution state by |0 so i) (i.e., ground 
state at t = 1). At the beginning of the algorithm, 
at time t — 0, the solution has approximately equal 
small overlap with all eigenstates, | (<?!>soi IV'-í (0)) | 2 ~ 
The solution probability is therefore distributed over all 
eigenstates. During the evolution this probability grad- 
ually gets "concentrated" in the ground state, so that 
after we have passed the minimal gap A at í m i n we have 
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K^soll^Mí > ímin)| 2 ~ 1 While |(0 so l|-0i>o(í > ímin)| 2 ~ 

0. With time the probability therefore "flows" towards 
the ground state. So even though the RMT core occupies 
high energies it could be important for the adiabatic al- 
gorithm because the solution probability for small times 
is found also in these high energy eigenstates. To check 
how much the RMT core states participate in this prob- 
ability flow we have calculated the total probability at 
time t to find the solution \4> S o\) in the eigenstates higher 
than ï-th, 

JV-l 

p( ï ,t)=J2\(^ ] (t))\ 2 - (7) 

Due to the normalization we of course have p(0,t) = 1, 
while at the end of the algorithm all probability is in the 
ground state, p(0, 1) = l,p(i > 0, 1) = 0, due to our defi- 
nition of the final ground state being the solution. Note 
that p(i, t) gives the failure probability if we stop the al- 
gorithm at time t and do not extend it to final 4=1. 
The failure probability at the end, at t = 1, provided 
we e.g., remove all levels higher than z-th from time i 
onwards is not necessarily equal to p{i,t). It is proba- 
bly correlated with p(i,t) but with details depending on 
how we "remove" the levels at time t. In Fig. ^| we 
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FIG. 10: (Color online) Logarithm (base 10) of the probability 
p(i, t) [Eq. J7J] to find the solution in eigenstates higher than 
i-th. Four black curves are isolines at 10 _1 , 10 -2 , 10 -3 , and 
10 -4 probability. The black ellipse has the same location and 
the size as in Fig. [5] and the 3-SAT instance is also the same. 



show the dependence of p(i,t) for the same test 3-SAT 
instance with n = 10 variables used before. We can see 
that the RMT core states are of limitcd importance. For 
instance, if we are satisfied with the error probability of 
10% (i.e., the probability to find the solution of 0.9), for 
times larger than t w 0.4 the states higher than i « 300 
(location of the isoline 10 _1 ) are not important as long as 
the probability flow for levels with i < 300 stays the same. 
Fig.^|nicely illustrates that as time progresses the high 
energy states become less and less important. This hap- 
pens already before the actual minimal gap is reachcd at 
ímin ~ 0.56. One consequence of this is that the adiabatic 
algorithm is expected to be more insensitive to the cou- 
pling of high energy eigenstates to the environment. For 
some results regarding the stability of quantum adiabatic 
algorithms to perturbations, see |29ll3Ct l3l|. 



VI. CONCLUSION 

We have shown that the properties of the ground state 
of the adiabatic algorithm for single-solution 3-SAT in- 
tances are very different from those of random vectors 
occurring in RMT. The distribution of gaps between the 
ground state and the first excited state for random 3-SAT 
problems with one solution shows a transition from GOE- 
like to Poisson-like distribution with increasing problem 
size. What is more, the distribution obtained for n = 18 
does not yet seem to be the asymptotic one. The ground 
state is also relatively weakly entangled as compared to 
RMT predictions. On the other hand, the central bulk 
portion of the spectrum is well described by RMT but 
has a limited influence on the flow of probability to the 
ground state. Therefore, RMT seems to be of limited use 
in describing Standard quantum adiabatic algorithm for 
single-solution 3-SAT instances. 
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